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Abstract 

The Poisson structures for 3D systems possessing one constant of 
motion can always be constructed from the solution of a linear PDE. 
When two constants of the motion are available the problem reduces to 
a quadrature and the structure functions include an arbitrary function 
of them. 


1 Introduction 


Three dimensional dynamical systems have deserved much attention both in 
view of their intrinsic mathematical relevance and of their wide interest in 
domains such as Mechanics |I], optics i-i. dynamic of interacting popu¬ 
lations ||-|[y modeling of fluid turbulence [§|-[|§, wave interaction models 
0-0- dynamo theory [|13|], and several other areas of physical, chemical 


or biological importance. A more recent related issue concerns the Poisson 
structures of these systems. This question has been contemplated both from 
the point of view of their existence [14] and of their explicit determination 


nj-p 


As a rule, 3D systems possess a Poisson structure whenever a suffi¬ 
cient number (e. g. two ) of independent constants of motion exist 


IS . Their 


explicit construction however is a partly open problem which, as we show 


1 





























in this letter, can always be solved when two constants of the motion are 
known. 


In a recent paper Gumral and Nutku [TO reduced the problem of the 


determination of the structure functions associated to the Poisson structure 
of 3D dynamical systems to the solution of a quasi-linear partial differential 
equation (eq. (70) of their paper) in three independent variables. When 
two independent constants of the motion are known, the problem therefore 
reduces to a Ricatti equation. In this letter we show that in fact the Poisson 
structures of 3D systems possessing one constant of the motion can always 
be obtained from the solutions of a linear partial differential equation and 
that when two independent constants of the motion are known the problem 
reduces to a quadrature. In such cases the resulting structure functions 
may involve an arbitrary function of the constants of the motion. When 
only one constant of the motion is available, the problem can frequently be 
handled in the sense that a generalized Hamiltonian formalism may still be 
constructed in terms of some particular solution of the pertaining equations. 
To illustrate the procedure we shall consider examples of both types. An 
interesting example with only one constant of motion is a five parameter 
version of the 3D Lotka-Volterra system for which a Hamiltonian formalism 
can be constructed. A four parameter version of the same system was already 
known to possess a bi-Hamiltonian structure H . 


2 Generalized Hamiltonian Structures 

In this section we consider a generic dynamical system in N dimensions 

= u M (x, t) , fM=l,...,N, (1) 

where is a sufficiently smooth vector held (in general E C°°), x = 
(x 1 ,^ 2 ,... ,x N ), and the over dot denotes derivative with respect to t. In 
addition we consider a function if(x, t) satisfying 

dH _ dH 
dt dt ’ 

along any phase trajectory, that is 

v> i d ll H = 0, 


( 2 ) 

(3) 
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where d fl indicates the partial derivative with respect to and repeated in¬ 
dices represent the Einstein’s summation convention. In specific applications 
the function H is typically a time-independent first integral of ([Ij) valid over 
a region in phase space. Finally we consider an anti-symmetric matrix J 
which satisfies the Jacobi identities HI 


= 0, (4) 

and therefore provides a generalized definition for the Poisson bracket 


[F, G\ = d^F d v G , 


(5) 


of functions F and G in phase space. 

Definition: System (|T|) is said to be Hamiltonian iff there exists a func¬ 
tion FI satisfying (|2|) and an anti-symmetric matrix J satisfying the Jacobi 
identities such that 

v 11 = ,P W d u H . (6) 

In such case, H is called the Hamiltonian of (:1) and J the associated Lie 
tensor or matrix of structure functions. 

For given H(x,t) let us label the variables such that d^H fy 0. 

Lemma: if for some anti-symmetric matrix J IW and //(x. t) satisfying ([|), 

v s = J™ d u H , for 8 = 1 ,... , (N - 1) (7) 

then 

v N = J Nv d v H , for s = N. (8) 

The implication of the lemma is clear: if a dynamical system possess a con¬ 
stant of motion and N — 1 of its equations are in a Hamiltonian-like form 
then its last equation is necessarily of the same form. 

Proof: Use ([/]) to expand (J3,) in terms of J sv 

N -1 

£ d s H J sv d u H + v N d N H = [H, H] + d N H (y N - J Nv d v H ) = 0. (9) 

•S=l 

The proposition now follows from the last equality, the property of the Pois¬ 
son bracket, and the fact that d^H fy 0. 
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Corollary: When ([7j) holds, (N — 1) components of the Lie tensor can be 
represented in terms of the remaining (N — 1 )(N — 2)/2 ones: 


j N /i _ _j/UV 


x \ 

f- E J 1 " d v H) / (d N H) 

v=l ) 


( 10 ) 


The proof is constructed by solving (0) for J Nv . 

Thus, any function satisfying (pi) recast system (HD in the “pre-Hamiltonian’' 
form (| 6 |). The algorithm consists of taking (N — 1 )(N — 2)/2 arbitrary func¬ 
tions J' w (// < v = 1,..., (N — 1)) and filling up the Lie matrix with them, 
their anti- symmetric counterpart and the J Nfl given by (|T0j). To complete 
the hamiltonization process we finally demand that J^ v obey the Jacobi iden¬ 
tities which then become their determining equations. 

In 3D the Jacobi identities are compatible with a generic conformal rescal¬ 
ing of the structure functions |T5|-|[T6j. Such invariance is, however, partially 


restricted for J^ u constrained by relation ([[(]). In fact, when a constant of 
the motion exists, the generic conformal invariance admitted by the Jacobi 
identities in 3D become a scale invariance by functions of the constants of 
motion and the Casimir only. We now prove this statement: 

Suppose that = 7 ( 2 ;) J 7 "' are structure functions for some Hamilto¬ 
nian H when are the corresponding structures functions for H. The 
equivalence of the two representations for the same system implies 


J l,v {~ld v H - d u H ) = 0 . 


( 11 ) 


As a Hamiltonian for a 3D autonomous system, H is necessarily a function of 
the two independent constants of the motion. Locally these constants 
may be chosen as the Hamiltonian H and the Casimir of -A". Therefore, 
without loss of generality, 


H = F(H, C). (12) 

Replacing H = F(H, C) in (|TTD and considering that the Casimir commutes 
with all functions of the dynamical variables yields 

7 = i/fl = C). (13) 

To close the section we remark that in general the Jacobi identities form 
an overdetermined system of non linear equations in the unknown . In 
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fact, when a constant of the motion exists they form a set of N\/ (3 \{N — 3)!) 
equations in (N — 1 )(N — 2)/2 unknowns { J\ fi < u — 1,, (N — 
1)}. For N < 3 the Jacobi identities are automatic and the hamiltonization 
process is trivial. For N > 3 the resulting system of equations is usually 
overdetermined. For N = 3 (which is the object of this letter) they become 
a linear partial differential equation. When two constants of motion are 
available, the problem is reduced to a first order ODE with solutions in terms 
of quadratures. The general solution for the structure functions contains an 
arbitrary function of the constants of the motion. The prove of this statement 
is presented in the following section. 


3 Three dimensional systems 

For N = 3, equations (0) read 

v 1 = J 12 d 2 H + J 13 d 3 H, (14) 

v 2 = -J 12 d 1 H+ J 23 d 3 H. (15) 

We can solve these equations in terms of some function J to be determined 
latter 


J 12 

J 13 

J 23 


J, 

(v 1 -Jd 2 H)/(d 3 H), 
(y 2 + Jd 1 H)/(d 3 H). 


(16) 


Recall that a different labeling of the variables must be adopted when d 3 H = 
0. This implies that equations (|T6D undergo a cyclic permutation of the 
indices (1 —> 2 —> 3 —> 1) that exchanges d 3 with either d 2 or d\. This 
permutation must be applied simultaneously to equations (|T8|-p~^) bellow. 


Equations (0) express all the components of the Lie tensor in 3D in terms 
of the symbol J determined by the Jacobi identity. To see this we substitute 
( |I 6 |) into (d) and multiply the result by d 3 H. The cross derivative and the 
quadratic terms cancel out, and, after using (|j) and some simple algebra, we 
obtain 

v^d^J = AJ+B, (17) 
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where 


(18) 


A = fyv" - 


(d 3 iA)(d^H) 

d,H 


and 


B = 


A d 3 v 2 — v 2 <9,3n 1 

dJ1 ' 


(19) 


In this equations the indices 1, 2 and 3 are linked to the indices in equa¬ 
tion (|l6| ) and, as already mentioned, should undergo a cyclic permutation 
whenever d :i II = 0. 

Equation (0, which is the Jacobi identity in terms of J, is the key 
equation in the solution of the hamiltonization problem of 3D systems. As 
mentioned before, this last condition is a linear Erst order PDE. When two 
constants of the motion are known, one of the characteristic equations for (|T7|) 
becomes a linear Erst order ODE (two out of the three dynamical variables 
are solved in terms of a third one and the constants of motion). In this case 
J can be solved by quadratures and the structures functions may always 
involve an arbitrary function of the two constants of the motion. 

It is important to remark here that the scale invariance referred to before 
is preserved by equations ((T7| - [T9|) : if H is a Hamiltonian with associated Lie 
tensor (JT(3|) (expressed in terms of J) then F(H,C), F arbitrary, is another 
valid Hamiltonian but with associated Lie tensor J = J)(dF/dH). In fact, it 
can be easily checked that if J satisfies equation (|T7j) for A and B calculated 
from LJ, then J satisfies the same equation with A and B calculated from 
F(H, C). The multiplicity of solutions for the hamiltonization process is still 
greater if we consider that equation (p~7|), being a first order linear PDE, 
in general admits an infinity of solutions. Sometimes, as we shall show in 
section 4, solutions in terms of arbitrary functions can also be found. 

As final remark of practical interest we point that any particular solu¬ 
tion to (jT7j) provides a non trivial Poisson structure for the system under 
consideration. Interesting particular cases are obtained when the system 
has B — 0, in which case J = 0 is the simplest possible solution or when 
H(x, t) = —f(H, C) A(x, t) in which case J = f(H,C ) is an equally simple 
particular solution. 


6 







4 Sample applications 

To illustrate some of the various possibilities offered by the procedure pro¬ 
posed above we consider now three sample problems in detail. Several other 
systems can be equally treated by repeating a similar sequence of steps. A 
more exhaustive and detailed list of examples is being prepared and will soon 
be submitted for publication elsewhere. 


4.1 The ice skate problem 


As a first example we consider the ice skate problem studied by Lucey [17 
In the original treatment one of the equations was x 2 = 0 which can be solved 
for x 2 = —a = constant. This transforms the 4D system into a 3D one. We 
next relabel the variables using the replacements x 4 —> x 3 —> x 2 . In this 
notation the ice skate system reads 


x 1 = 


x 2 = x 3 


x 3 = ax 3 tanx 1 . 


( 20 ) 


We now verify that H\ = x 3 sec(x 1 ) is a constant of the motion and therefore 
a candidate as the Hamiltonian for the system. Use of this Hamiltonian as a 
source in (p~8| ~ [T^ ) gives A = 0 and B = —acosx 1 . We next verify that one of 
the characteristic equations of (|T7|), 

dJ 




dx 1 


= cosx 


is separate from the others and can be integrated in the form 

J = smx 1 + F{H 1 ,H 2 ), 


( 21 ) 


( 22 ) 


where F is an arbitrary function of constants of the motion. The following 
structures functions are now obtained by substitution of ([221) hi 


ju = _j 2i = sinx i + ^ 

J 13 = — J 31 = — a cosx 1 , (23) 


J 23 = —J 32 = x 3 (secx 1 + F(Hi, H 2 ) tanx 1 ) . 

To find a Casimir of the algebra it is necessary to specify the function F. In 
particular for F — 0 , C\ — ax 2 + x 3 tanx 1 , is a Casimir of the system. 
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We now verify that the ice skate problem is in fact completely integrable 
since C\ is another constant of motion, functionally independent of H\. We 
may therefore use H 2 = C\ as an alternative Hamiltonian in which case 
A = B = —a cot x 1 and 


J=-l + F{H 1 ,H 2 ) S mx 1 , (24) 

for any arbitrary function F of the constants of the motion. This generates 
the alternative Poisson structure 

J 12 = _J 21 = _! + f(H 1 ,H 2 )sinx\ 

ji3 = _J3i = ,H 2 )cosx\ (25) 

J 23 = — J 32 = x 3 f«- tan a ; 1 + F(Hi, H 2 ) sec a: 1 ) . 

For F — 0 the new algebra has Casimir C 2 — H\ — a ; 3 sec ad. It is worthwhile 
remarking that the algebras corresponding to J and J are independent and 
not connected by conformal transformations. The rescaling invariance by 
functions of the Hamiltonian and the Casimir can still be applied to generate 
families of equivalent Poisson structures. 

We have presented a complete and novel solution to the hamiltonization 
of the ice skate system. This problem illustrates clearly the basic features 
of the routine proposed. In the next subsection we analyze another sample 
problem and show that a well known system may still exhibits some novel 
and perhaps unexpected features. 

4.2 The Euler top 

We now revisit the Euler top problem and show that also this system pos¬ 
sesses families of alternating Hamiltonian and Casimir with the associated 
Poisson structure involving arbitrary functions of the two. 

The classical Euler top system |]J can be written as 

x 1 = {h - h)x 2 x 3 / (J 2 / 3 ), 
x 2 = {h - Ii)x 3 x 1 /(I 3 Ii) , 
x 3 = {I\ - I 2 )x 1 x 2 /(I 1 I 2 ) , 

8 


( 26 ) 


(27) 


and admit the (kinetic) energy 


H 


1 

2 


^•’) 2 


+ 




and the angular momentum 


l = (x 1 ) 2 + (x*y + (x 3 y, 


(28) 


as independent constants of the motion. Taking the energy as the Hamilto¬ 
nian we find B — 0 and 


A — (Ii — I 2 )x 1 x 2 /( I\I 2 x 3 ) = x 3 /x 3 . (29) 

For these values of A and B the solutions J to the basic equation (|T7|) are 
easily found and the corresponding structure functions become 

J 12 = -J 21 = -x 3 (l + F(H,L)) , 

J13 = _J31 = x 2 (1 + I 3 F(H, L)/I 2 ) , (30) 

J23 = -J 3 * = -x 1 (l + I 3 F(H,L)/I 1 ) , 

where again F is an arbitrary function of its arguments. When F — 0 the 
usual results are recovered with the angular momentum L as the Casimir. 

Alternatively, starting with the angular momentum L as the Hamiltonian, 
leads to the same forms for A and B but the new structure functions become 
(where F is also an arbitrary function of its arguments) 

ji2 = -jn = x z {\/h + F{H,L)) , 

J 13 = - J 31 = -x 2 (l // 2 + F(H, L)) , (31) 

J 23 = - J 32 = x 1 (l/h + F(H, L)) . 

For F = 0 the kinetic energy H is the Casimir of the new algebra. 

The analysis presented here shows that the Euler top system admits a 
continuous family of Poisson structures in which the kinetic energy and the 
angular momentum play alternatively the roles of Hamiltonian and Casimir. 
Each pair of Hamiltonian and Casimir have associated Lie algebra that are 
independent, that is, that are not connected by conformal rescaling. 
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4.3 The 3D Lotka-Volterra system 

The 3D Lotka-Volterra system and some of its special subsystems play an 
important role in modeling many physical, chemical and biological processes. 
Their associated vector field is defined (in their most general form with Ver- 
hulst terms b„ y 0 ) by 

v k = x k (a k + b^) k = l,...N. (32) 

In equation (|32]) and throughout the rest of this letter sum over the Latin 
index k is not implied. 

Cairo and Feix Q found several invariants or first integrals for N-dimensional 
Lotka-Volterra systems. In particular for N = 3 and det(fejj) = 0, one such 
invariant is 

H = H{x\ x\ x 3 ) e~ st = {x^ixY^ye-^ (33) 

where a, /3, 7 and s are given by 

a = &22&31 — b 2 ib 3 2 , 

P = biib 32 ~ bi 2 b 3 i , 

7 = & 12&21 — ^ 11^22 j 

s = Oja + a 2 P + 037 . 

Under the constraint s = 0, that is, for 

01 (^ 22^31 — b 2 ib 32 ) + 02 (^ 11^32 — b\ 2 b 3 \) + 03 (^ 12^21 — bnb 22 ) = 0 , (34) 

the function H becomes a time-independent first integral. Bellow we present 
the derivation of a Poisson structure for this system under the less restrictive 
conditions that we can find. This will ultimately imply that among the ini¬ 
tially free coefficients and bij only five remain arbitrary. Unfortunately no 
other constant of the motion is known for the same number of free parameters 
and complete integrability cannot be sought. It is however interesting to ver¬ 
ify that under appropriate conditions one can construct a Poisson structure 
without being actually able to integrate the equations completely. 

In order to tackle the basic equation ([L7]) we first calculate A and B in 
terms of the vector field v M and of the constant of the motion H. For this, 
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(35) 


equations (|T8| - |T9D and ( |33| ) yield 

A = ai + a 2 + bnx 1 + 6 22 x 2 + 6 33 x 3 + ( b lfl + 6 2 /( )x /i , 


B — y-^[(ai&23 — a 2 b\ 3 ) + (b 23 bi^ — bi 3 b 2 ^)x^] , 


(36) 


where we introduced the symbol U = x x x 2 x 3 to simplify the notation. To 
proceed we write the characteristic equation associated to 0) 


dx k 


dJ 


x k (a k + b kfl x dJ + B’ 


which after some element but tedious algebra implies 

dU 


dJ 


U [A + a 3 + ( 6 3 i — 6n)x 1 + (632 — 6 2 2)x 2 ] did 4 - B 


(37) 


(38) 


We now multiply the numerator and the denominator of the first term in 
equation (|38D by e/^H where e is an arbitrary constant to be determined 
later, and re-write the characteristic equations in the form 


d x k 


dlJ-eU/^H)} 


x“(a t + h^) A[J — eU/^H)] - B' (U/(jH)) ’ 
with B' dehned by 


(39) 


B' = 03623 - a 2 b l3 - ea 3 + [e(b n -b 31 ) + 623611 - 613621 ] x 1 

+ [ 4(622 632 ) 3-623612 613622 ] x 2 . (40) 


Equation ([39]) is difficult to treat in its general form. One particular solution 
however can be readily found by imposing the additional condition B' = 0. 
As can be easily checked, under this condition 

J=^r = -(x 1 ) 1 -^ 2 ) 1 -^ 3 ) 1 - 7 - (41) 

7 H 7 

satisfies equation ([59]) and therefore the fundamental equation ([L7]). The 
condition on B' produces the value of the arbitrary constant e and (for arbi¬ 
trary x 1 and x 2 ) two additional constraints on the coefficients of the system. 
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This represents a total of four constraints (det( 6 y) = 0, s = 0 and two out 
of conditions (42- 44]) below) on the twelve initially free parameters. This 
corresponds to a total of eight free parameters which in fact reduces to five 
since three parameters can always be eliminated by rescaling. 

Condition (|4(J) for arbitrary x 1 and x 2 , implies three equations one of 
which can be solved for e when either (i) a 3 7 ^ 0 or (ii) 6 31 7 ^ bn or (in) 
632 7 ^ & 22 - These conditions imply the following relations 


(cq&23 — o, 2 bi 3 )(b 3 i — &11) = ^3(^23^11 — ^1.3^21), (42) 

( 01&23 ^ 2 ^ 13 ) (^32 b 22 ) = o, 3 (b 23 bi 2 — b± 3 b 22 ) , (43) 

(^23^11 — ^13^21) (^32 — ^22) = (^23^12 — ^13^22) (^31 — ^11) , (44) 


which we interpret as follow: 

i) when a 3 7 ^ 0 apply (fI2| - [l5| ) and use 

e = -(aib-23 ~ a 2 b 13 ) ; (45) 

03 

and/or 

ii) when b 31 7 ^ bn apply i |T 2 | . [Ti] j and use 

£ = (&23&11 — &13&2l))/(&31 — Wi ); (46) 

and finally 

iii) when b 32 7 ^ b 22 apply (|43| - |4^) and calculate e from 

e = (&23^12 — ^13^22)/(^32 _ b 22 ) . (47) 

To complete the calculation insert e in (|4l|) and substitute J into (|T6|). 
This provides a Poisson structure for the 3D Lotka-Volterra system, namely: 

ji2 = _j2i = J_i x 2 x * (48) 

7/1 

^ (y - ^) . ( 49 ) 

j» = -j« = ?*(y + fEL x >). (50) 

In Nutku’s analysis H of the Lotka-Volterra system, one of the first inte¬ 
grals was the logarithm of the function used above and no Verhulst (diagonal) 
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terms were included. Like in the present treatment the determinant of the 
coefficient of the quadratic terms was set to zero and an additional constraint 
was imposed similarly to condition ( |33| ) . This more constrained system (four 
free parameters) has a known second constant of motion (or a constant of 
motion and a Casimir) and consequently admits a bi-Hamiltonian structure. 
The component J 12 of the Lie tensor in Nutku’s analysis satisfies condition 
( pT|) but is a different particular solution. 


5 Conclusions 


We have presented a procedure for constructing the Poisson structure of 
dynamical systems possessing a constant of motion. In three dimensions 
the problem is reduced to the solution of a linear PDE. When two time- 
independent first integrals of the system are known the procedure can be 
applied twice and yields bi-Hamiltonian structures involving arbitrary func¬ 
tions of the constants of motion (or a constant of the motion and a Casimir). 
The technique was applied to a few sample systems to show how it oper¬ 
ates. In particular we showed that the ice skate system analyzed by Lucey 
IDU is completely integrable, that the Euler top systems admits an infinite 


family of bi-Hamiltonian structures and that a five free parameter version 
of the 3D Lotka-Volterra system can be cast in a generalized Hamiltonian 
form. Other 3D systems for which one or more time-independent constants 
of motion are known can be equally treated. This is being done for several 
rescaled version of such systems found in the literature. The results are still 
in preparation (some preliminary results can be found in ref. |18|) and will 
be soon submitted for publication elsewhere. 
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